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The dynamic behavior of polyelectrolyte chains in the oligomer range is investigated with coarse- 
grained molecular dynamics simulation and compared to data obtained by two different experimental 
methods, namely capillary electrophoresis and electrophoresis NMR. We find excellent agreement of 
experiments and simulations when hydrodynamic interactions are accounted for in the simulations. 
We show that the electrophoretic mobility exhibits a maximum in the oligomer range and for the 
first time illustrate that this maximum is due to the hydrodynamical shielding between the chain 
monomers. Our findings demonstrate convincingly that it is possible to model dynamic behavior 
of polyelectrolytes using coarse grained models for both, the polyelectrolyte chains and the solvent 
induced hydrodynamic interactions. 
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INTRODUCTION 

Electrophoresis methods are widely used to separate 
(macro)biomolecules P,Q such as peptides, proteins, and 
DNA, as well as synthetic polymers @, 0]- Short poly- 
electrolytes (PEs) can conveniently be separated in free 
solution without the aid of a gel by capillary electrophore- 
sis (CE). Additionally, CE is employed to characterize 
the hydrodynamic properties of charged biomolecules, in 
particular the electrophoretic mobility, u, the diffusion 
coefficient, and the hydrodynamic radius @,@j. 

Alternatively, these properties can be determined by 
pulsed field gradient (PFG) NMR @, @]. With a com- 
bination of diffusion NMR and electrophoretic NMR, 
the charge of macromolecules the influence of ionic 
strength and the dielectric constant of the solution have 
been investigated (loj . 

Several studies of PEs of well defined length in the 
short chain regime have shown that the free-solution mo- 
bility /i exhibits a characteristic behavior 13, 
from the monomer mobility, /j,q, onwards, /i increases 
towards a maximum, /i max , that occurs for chains of a 
specific degree of polymerization, N max . After this maxi- 
mum, /it decreases slightly to reach a constant value, fjtpo, 
the so-called free draining mobility. Very little is known 
about the origin of the maximum, yet the knowledge of 
the precise dynamical behavior of a PE is a prerequisite 
for designing possible applications in microfluidic devices, 
such as electrophoretic separation or gene sequencing. 

Whether or not a mobility maximum appears and the 
precise value of -/V max seems to depend on the flexibil- 
ity of the chain. While a maximum has been observed 



for flexible PEs such as single stranded DNA or sul- 
fonated polyst yre ne (PSS) at values around 10 repeat 
units for A^ayjTll, [tH. [l3j. semi-flexible double-stranded 
DNA does not exhibit it. The small number for 7V max and 
the strong influence of the persistence length shows that 
a short chain behavior must be responsible for this. Fur- 
thermore, it has been shown that the /i max is diminished 
under elevated salt conditions. 

Existing theories 11 , ID, 14 , li|, lj| have been success- 
ful in describing the qualitative behavior of an initially 
rising mobility as well as the constant long-chain limit, 
but they have not been able to reproduce this maximum 
or explain its origin. 

To provide a fundamental understanding of the in- 
volved dynamics of PEs, we propose a coarse-grained 
Molecular Dynamics model that takes all charged par- 
ticles (i.e. PE repeat units, counter ions, and additional 
salt) and electrostatic and hydrodynamic interactions 
(HI) between them into account. In the following we 
demonstrate that this model quantitatively reproduces 
experimental results obtained by two completely differ- 
ent experimental techniques. In addition we suggest a 
microscopic interpretation of the size-dependence of fi 
for short PE chains, based on hydrodynamic shielding, 
which gives fundamental insight into the interplay of hy- 
drodynamic friction and charge correlations for charged 
macromolecules . 



SIMULATIONS 

Model. We simulate a flexible PE using a bead-spring 
model employing the ESPResSo package All pa- 
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rameter values are chosen to match the properties of PSS 
used in the experiments and given in reduced units with 
kT = 1.0 and a = 2. 5 A being the energy and the rel- 
evant length scale. The beads (chain monomers) are 
connected by FENE bonds with stiffness k = 30, and 
maximum extension R — 1.5. Additionally a truncated 
Lennard-Jones or WCA potential with depth e = 0.25 
and width a = 1, is used for excluded volume interac- 
tions. Each monomer has a charge of q = — 1 in units 
of e. Monovalent counter ions, q — 1, and monovalent 
positive and negative salt ions are subjected to the same 
LJ potential giving all particles the same size. The sim- 
ulations are carried out under periodic boundary condi- 
tions in a rectangular simulation box of size L — 34 (for 
N = 2) to L = 89 for (N = 32) resulting in constant 
monomer concentration of approximately 1 g/1 or 5 mM. 
All electrostatic interactions are calculated with the P3M 
algorithm. The Bjerrum length l B = e 2 /47re/sT = 2.84 
in simulation units corresponds to 7.1 A (the Bjerrum 
length for water at room temperature). Together with 
the model's average bond length of b = 0.91, we com- 
pute a Manning factor of£ = ?s/& = 3.12. The inclusion 
of HI is done via frictional coupling of the beads to a 
Lattice Boltzmann (LB) fluid as detailed in [lij . The 
modeled fluid has a kinematic viscosity v = 3.0, a fluid 
density p = 0.864, and is discretized by a grid with spac- 
ing a = 1.0. The coupling parameter is Tbarc = 20.0. 
The simulation time step is r = 0.01. 

To determine the impact of HI, we compare the results 
to simulations with a simple Langevin thermostat that 
does not recover long-range hydrodynamic interactions 
between the monomers, but only offers local interaction 
with the solvent. 

Transport coefficients. We determine two differ- 
ent transport coefficients for the model PE that are like- 
wise determined in the associated experiments. The sin- 
gle chain diffusion coefficient is obtained from the cen- 
ter of mass velocity auto correlation function: D = 

In CE experiments, the electrophoretic mobility [i of 
the solute is determined by /i = v/E = Ll/Vt, where v 
is the velocity, E is the electric field, V is the applied 
voltage, L is the total length of the capillary, I is the 
migration (or effective) length up to the detector and t 
is the detection time of the solute. 

In the simulations, we use a Green-Kubo relation to 
obtain fi at zero electric field. This approach has been 
successfully applied in simulations to determine the elec- 
trophoretic mobility of charged colloids ([l3|). The chain 
mobility is calculated from the correlation function be- 
tween the center of mass velocity of the polyelectrolyte 
chain and the velocities of all charged particles in the sys- 
tem: fi = j^J^tli fo°(fii(fy ' v c (t))cLt. This method 
guarantees that no conformational changes of the chain 
structure or the ion distribution are induced by an artifi- 
cial high external field, which is sometimes used in other 



simulations to separate the directed electrophoretic mo- 
tion from Brownian motion within reasonable computing 
time. Our method enables us furthermore to obtain both 
transport properties from the same simulation trajecto- 
ries without additional computational effort. 21 



EXPERIMENTS 

Capillary electrophoresis (CE) is an analytical 
separation technique based on the differential migration 
of ionic species under electric field [H . 

The CE experiments were performed using an Agi- 
lent technologies capillary electrophoresis system (Agi- 
lent, Waldbronn, Germany). Capillaries of 33.5 cm (25 
cm to the detector) length, and 50 /im diameter were 
prepared from bare silica tubing purchased from Supelco 
(Bellefonte, PA, USA). New capillaries were conditioned 
by performing the following flushes: 1M NaOH for 30 
min, 0.1 M NaOH for 10 min, and water for 5 min. 
Samples were introduced hydrodynamically (~ 4 nL) at 
0.5 g/L concentration (~ 2.5 mM monomer concentra- 
tion). The electrolyte was pure water. Solutes were de- 
tected at 225 nm. The electric field was kept constant 
at 224 V/cm (V=+7.5 kV). The polarity of the applied 
voltage on the inlet side of the capillaries was positive. 
All the experiments were performed at 27 °C. Electro 
osmotic mobilities were determined from the migration 
time of a neutral marker (mesityl oxide, ~ 0.1% (v/v) 
in the electrolyte). The sodium polystyrene sulfonate 
standards (M^l.430 x 10 3 , 5 x 10 3 , 8 x 10 3 , 16 x 10 3 , 31 x 
10 3 ,41 x 10 3 ,88 x 10 3 ,177 x 10 3 ,350 x 10 3 ; M w /M n ca 
1.1) were purchased from American Polymer Standards 
Corp. (Mentor, OH, USA). These standards have ver- 
ified sulfonation rates larger than 88% and are almost 
fully charged. Purified water delivered by an Alpha-Q 
system (Millipore, Molshcim, France) was used to pre- 
pare all electrolytes and sample solutions. 

In electrophoretic NMR, diffusion and elec- 
trophoretic motion are separated by the design of the ex- 
periment [2(j. No gel has been used, so that self diffusion 
and free electrophoresis are measured. The PSS sam- 
ples have been obtained from Fluka. To minimize the ef- 
fects of variations of the ionic strength [l(| , samples have 
been dialyzed against water (cut off volume 0.5kDA), and 
subsequently dried under vacuum. For all experiments 
a monomer concentration of 5 mM in deuterated water 
has been used. The diffusion experiments have been per- 
formed on a Bruker Avance 500 NMR spectrometer op- 
erating at a Larmor frequency of 500MHz for protons 
equipped with a DIFF30 probe head generating a maxi- 
mum pulsed field gradient strength of 12T/m. The gradi- 
ent pulse duration and diffusion times have been adjusted 
between 0.8 to 2ms and 8 to 20ms respectively for optimal 
resolution for each molecular weight resulting in different 
diffusion coefficients. Because of their narrow molecular 
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weight distribution, diffusion coefficients have been deter- 
mined by a linear fit to the Stejskal- Tanner equation Q. 
Electrophoresis NMR experiments have been performed 
on a Bruker Avance 300 NMR spectrometer operating at 
a Larmor frequency of 300MHz for protons with an in- 
house-built electrophoresis probe head utilizing a Bruker 
micro2.5 imaging gradient system generating magnetic 
field gradient strength of up to lT/m. The flow times 
varied between 15 and 50 ms and the gradient pulse du- 
ration between 3 and 6 ms respectively with gradient am- 
plitudes between 0.3 and 0.6 T/m. The electric field has 
been linearly incremented between -140 V/cm and +140 
V/cm. The electrophoretic mobility is a model- free read 
out from the two-dimensional electrophoresis NMR spec- 
trum correlating the chemical shift, that identifies the 
moving species, with electrophoretic mobility Q. 
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FIG. 1: The normalized diffusion coefficient, D/Do, for PSS 
of different lengths N. The diffusion for simulations with HI 
and the experimental data agree and show a scaling exponent 
of m — 0.63. Without HI the diffusion scales with m = 1. 



Diffusion. Recently, Stellwagen et al. [13j observed 
that the diffusion coefficient, D, of PSS can be approx- 
imated by a power law scaling D = D N m , where Dq 
is the monomer diffusion coefficient, and m is the scal- 
ing exponent. In Figure [T] we compare the diffusion co- 
efficient obtained from simulations to the results from 
the NMR study. The simulated data is normalized by 
Do = 0.052 as obtained from the power law fit, and the 
experimental data by the monomer diffusion coefficient 
of D = 5.710~ 10 m 2 /s. The simulations with HI re- 
sult in a scaling exponent of m = —0.63 ± 0.01, which 
is in good agreement with value obtained from experi- 
ments, to = —0.64 |ld| . and a previously reported result, 
to = —0.617, by [l3j. Only for the very short chains, 
(N < 5), is a deviation from the prediction observed and 



a higher diffusion coefficient found in the simulations. For 
intermediate chain length, the coarse-grained simulation 
model with HI, reproduces the experimentally observed 
behavior. 

Without HI, the chains show the expected Rouse dif- 
fusion with an exponent of to = 1.02 ± 0.02. This simple 
model is clearly not applicable to mimic the experimental 
behavior of short PE chains. 
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FIG. 2: Normalized electrophoretic mobility /i//ifd as a func- 
tion of the number of repeat units N for simulation with HI, 
from capillary electrophoresis (CE), and from electrophoretic 
NMR. The inset compares simulations with and without HI. 



Electrophoretic mobility (//). The results for mea- 
surements of /j, in pure water without additional salt are 
shown in Fig. [2] To account for the different viscosity of 
the solvents, we rescale the mobilities by the free-draining 
mobility /ifd as obtained for long chains. 

The experimental data sets agree within the accuracies 
of the individual methods and show the characteristic be- 
havior of the mobility in dependence to chain length. A 
mobility maximum for N = 10 is observed with capillary 
electrophoresis. This maximum for intermediate chains 
as well as the long chain behavior is successfully repro- 
duced in simulations with HI. For the first few oligomers, 
we observe a small difference which is in line with the 
deviation for the diffusion. 

On the other hand, as illustrated in the inset of Fig. [2l 
the simulation without hydrodynamic interactions fails 
completely to describe the short chain behavior and can 
only be mapped to the experimental data in the long 
chain limit. Therefore, we infer that the mobility maxi- 
mum can only be explained when taking into account HI 
between the PE and the surrounding solvent ,22] . 

To further illustrate this, we estimate the effective hy- 
drodynamic friction, r o fr, and the effective charge, Q e ff, 
in dependence of its length. During electrophoresis, the 
electrical driving force Fe = Q c sE is balanced by the 
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FIG. 3: Effective charge per monomer Q c s/N (squares) as 
estimated from counterion condensation and effective friction 
per monomer F e f{/N. Q e g is approaching Q/N = l/£ = 0.32 
as predicted by Manning theory. 



frictional force with the solvent Fp = T c sv, where v is 
the migration velocity induced by the electric field E. 
For the mobility, we then obtain: /i = Q c s/Y c s- The 
PE is surrounded by oppositely charged counter ions, 
some of which moving with the chain and thus reduc- 
ing its effective charge. We estimate this charge reduc- 
tion by subtracting the number of counter ions that are 
found within 2 units of the chain from the bare charge 
N . This estimate is used together with the obtained mo- 
bility to calculate the effective friction as shown in Fig. 
[3J The effective charge for long chains is in agreement 
with the Manning prediction at infinite dilution, yielding 
Q/N = l/£ = 0.32. From Fig. [3J we see the impact 
of HI that results in the mobility maximum. Additional 
monomers are partly shielded from the flow by the other 
monomers decreasing the effective friction per monomer 
with chain length. This shielding is effective on short 
length scales and leads to a stronger initial decrease of 
the friction than the counter ion condensation reduces 
the effective charge. This leads to an increasing mobil- 
ity for intermediate chains. While the effective friction 
levels off quickly, the effective charge per monomer keeps 
decreasing slowly, reducing the mobility and causing the 
maximum. For long chains, effective charge and effective 
friction per monomer become constant, leading to the 
well-known free-draining behavior. 



CONCLUSION 

We performed a detailed study investigating the dy- 
namic behavior of short PE chains via MD simulations 
of a coarse grained model and via two different experi- 
mental approaches. The results of experiments and sim- 



ulations can be quantitatively matched and agree with 
the existing theory and predictions, as long as the simu- 
lation model correctly includes long-range HI. A simula- 
tion model that neglects HI fails to reproduce the short- 
length scale behavior of the diffusion coefficient and of 
the electrophoretic mobility. 

To our knowledge we demonstrated for the first time, 
that the transport coefficient of short PEs can be quanti- 
tatively modeled by coarse grained simulation techniques. 
No chemical details are needed to explain the experimen- 
tal results. Our model allows to simulate time scales oth- 
erwise out of reach for atomistic simulations and provides 
a microscopic understanding of the observed maximum 
in fi. From our observation we conclude, that hydrody- 
namic shielding between the chain monomers is the major 
cause for this previously little understood behavior, and 
that it is vitally important to include HI to model such 
systems. 

Having a simulation model at hand, that confirms the 
experimental data, opens new possibilities of investigat- 
ing the electrophoretic behavior of short PEs, which so 
far has not been fully explained by the existing theories. 

We thank B. Diinweg, U. Schiller, and G. Slater for 
helpful remarks. Funds from the the Volkswagen founda- 
tion, the DAAD, and DFG under the TR6 are gratefully 
acknowledged. 
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